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Abstract 

We show that the smooth hybrid inflation is naturally realized in a framework 
of supersymmetric axion model. Identifying the Peccei-Quinn scalar fields as a part 
of the infaton sector, successful inflation takes place reproducing the amplitude and 
i— H spectral index of the curvature perturbation observed by WMAP. A relatively large 

axion isocurvature perturbation and its non-Gaussianity are predicted in our model. 
The saxion coherent oscillation has a large amplitude and dominates the Universe. 
. ^ The subsequent decay of the saxion produces huge amount of entropy, which dilutes 

^ unwanted relics. Winos, the lightest supersymmetric particles in this scenario, are 

produced non-thermally in the decay and account for dark matter. 
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1 Introduction 

The Standard Model (SM) in particle physics is confirmed after the recent discovery of 
the Higgs boson at the LHC [1]. However, the SM is not believed to be a fundamental 
theory because of some theoretical problems. One of these is known as the strong CP 
problem. Although quantum chromodynamics (QCD) allows the existence of the CP- 
violating term in the Lagrangian, the experimental test of the neutron electric dipole 
moment suggests that CP must be preserved with very high accuracy implying that the 
CP- violating term must be extremely small [2]. The only plausible solution to the the 
strong CP problem so far is the Peccei-Quinn (PQ) mechanism [3]. In the PQ mechanism, 
a global symmetry written as C/(1)pq is spontaneously broken and a Pseudo-Nambu- 
Goldstone boson arises accordingly : it is called axion. The axion is assumed to acquire 
the potential dominantly through the QCD instanton effect, which drives the axion to the 
CP-preserving minimum [4] . 

On the other hand, the recent cosmological observations, including the Wilkinson 
Microwave Anisotropy Probe (WMAP) observations [5] strongly support inflation in the 
very early Universe. Generally, inflation is considered to be driven by some scalar field 
called inflaton with almost flat potential. Such a scalar field is not predicted in the 
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framework of SM. It is also revealed that our present Universe is filled with the non- 
baryonic cold dark matter (CDM), whose density parameter is found to be [5] 



where h is the dimensionless Hubble parameter in units of H = 100/i km/sec/Mpc 
with Hq being the present Hubble parameter. The SM does not contain any candidate 
of the CDM, and hence we are forced to go beyond the SM. Supersymmetric(SUSY) 
theories generally predict many scalar fields with flat potentials one of which may be an 
inflaton, and SUSY protects flatness of the potential against large radiative corrections. 
Furthermore, the lightest SUSY particle is a good candidate for dark matter due to R- 
parity conservation. 

Motivated by these observations, we consider a SUSY axion model in order to account 
for inflation and DM with solving the strong CP problem in a unified framework. In the 
previous papers [7, 8], we studied the inflation and the following cosmological history in 
a SUSY axion model. We pointed out that the SUSY hybrid inflation [9] naturally takes 
place by identifying the PQ scalar fields as waterfall fields. To reproduce the observed 
curvature perturbation, the PQ scale must be of order of f a ~ 10 15 GeV. Generally, 
such a large PQ scale leads to the overproduction of axions because the amplitude of 
the axion coherent oscillation scales as the PQ scale. However, considering the post- 
inflationary dynamics, the saxion (the scalar partner of the axion) dominates the Universe 
and the subsequent decay of the saxion produces a huge amount of entropy. The axion 
density is diluted and, as a result, it can have a right dark matter abundance [10, 11]. 

In this paper, we consider a variant inflation model: smooth hybrid inflation model [12, 
13] in the framework of the SUSY axion model by modifying the PQ sector appropriately. 
The idea was already noted shortly in the previous paper [8] , but complete analyses were 
not yet performed. We will see that the successful smooth hybrid inflation takes place in 
a SUSY axion model. Furthermore, considering the post-inflationary dynamics, we show 
that the saxion coherent oscillation inevitably dominates the Universe and the subsequent 
decay of the saxion can produce a correct amount of dark matter. 

This paper is organized as follows. In Sec. 2, we introduce a SUSY axion model we have 
in mind. In Sec. 3, the smooth hybrid inflation model is reviewed. In Sec. 4, reheating 
after inflation and the dynamics of saxion and its cosmological implications are discussed. 
In Sec. 5, we show the constraint from the isocurvature perturbation and calculate the 
non-Gaussianity. We conclude in Sec. 6. 

2 A supersymmetric axion model 
2.1 The potential of the SUSY axion model 

In this section, we briefly review the SUSY axion model. We consider the SUSY axion 
model whose superpotential is given by 



fi CDM /i 2 = 0.1126 ±0.0036, 
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Table 1: Charge assignments on the field content. 

where S is a gauge-singlet chiral superfield and $ and ^ are the chiral PQ superfields 
which are also gauge- singlet. X and X are the chiral superfields interacting with the 
PQ field at tree level and they also interact with the minimal supersymmetric standard 
model (MSSM) fields through the gauge interaction, fi and M are some mass scales, A 
is a dimensionless coupling constant and n is an integer larger than or equal to 2. This 
superpotential has the global C/(1)pq symmetry and the global U(1)r symmetry as well 
as the discrete Z n symmetry. Charge assignments of respective fields are summarized 
in table 1. The PQ superfields contain the axion (a), saxion (a, the scalar partner 
of the axion) and axino (a, the fermionic superpartner of the axion). There are two 
representative axion models : one is the Kim-Shifman-Vainshtein-Zakharov (KSVZ) model 
(also known as the hadronic axion model) [14] in which X and X are additional heavy 
quarks, denoted by Q and Q. The other is the Dine-Fischler-Srednicki-Zhinitsky (DFSZ) 
model [15] in which X and X are identified as Higgs fields, H u and H^. 
In the global SUSY limit, the F-term scalar potential is calculated as 

2 ra 2 |S| 2 |M 2 ( n_1 ) 



„2 | W 
* ^ ][{2(n-l) 



X[i ,,-n iW + m, (3) 

where the scalar fields are denoted by same symbols as corresponding superfields. Here 
and hereafter, X and X are set to be zero due to Hubble- induced masses or thermal 
corrections because X = X = is en enhanced symmetry point. The global minimum of 
the potential is placed at 

S = and = (4) 
where the PQ breaking scale f a is given by 

f a =(liM n - 1 ) 1/n . (5) 

This indicates the existence of the flat direction along which the PQ fields do not feel the 
potential, ensured by the fact that the £7(1)pq symmetry is extended to the complex U(l) 
due to the holomorphy of the superpotential [16]. 

The flat direction is lifted by the low-energy SUSY-breaking mass terms 

Vsoft = c im l /2 \y\ 2 + c 2 m 2 3/2 \ty\ 2 , (6) 

where m 3 / 2 is the gravitino mass and c\ and c 2 are real, positive and 0(1) numerical 
constants. This stabilizes the radial component of the PQ fields, |\&| and \^\ at 

n 1/4 / \ 1/4 
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respectively. The saxion is defined as the deviation from the minimum along the flat 
direction. Near the minimum of the potential (7), the axion and the saxon is related to 
the PQ fields as 

where F a is defined as F a = \/v 2 + v 2 . The PQ fields obtain vacuum expectation val- 
ues (VEVs) and, since the C/(1)pq is anomalous under the QCD, the axion obtains an 
instanton-induced potential and solves the strong CP problem via the PQ mechanism. 

2.2 The decay of the saxion 

Because the decay of the saxion is an important ingrediate in the following discussion, 
we here summarize the decay rate of the saxion. The saxion has an interaction with the 
axion through the kinetic terms of the PQ scalar fields as 

|^*| 2 + |^| 2 = (l + ^a) {^-{d,af + \{d,a)^ + . . . , (9) 

where £ is defined as £ = (v 2 — v 2 )/F 2 which is generally of order unity unless v ~ v i.e. 
C\ ~ C2 [17, 18]. The decay rate of the saxion into the axion pair is derived as 

- e m ° (10) 



64tt F 2 



where m a is the mass of the saxion which is same order of the gravitino mass. 

In the KSVZ axion model, the main decay mode of the saxion into the MSSM particles 
is that into two gluons whose decay rate is calculated as 

2 3 

r - " s 171(7 

^ gg 327r 3 F a 2 ' 1 } 

where a s is the QCD gauge coupling constant. Another efficient decay mode of the saxion 
is that into two gluinos if it is allowed kinematically, whose decay rate is given by 

- J^rf , (12) 

where d is an 0(1) numerical constant for m CT ~ m 3 /2 [19]. In the KSVZ axion model, 
the decay mode into the MSSM particles may be subdominant and the dominant decay 
mode may be the two axion decay unless c\ ~ c 2 . Assuming the decay into two axion is 
suppressed, the decay temperature is calculated as 

m n \ 3/2 /10 15 GeV\ 



T^-lOOMeV . (13) 

1 lOOTeV/ V F a J 
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In the DFSZ axion model, there exists the tree level coupling of the saxion with the 
standard model Higgses, so the saxion decays dominantly into the Higgses. The decay 
rate of the saxion into the lightest Higgs boson pair is derived as 



a—thh 



16tt V m, 



(14) 



where \x is the higgsino mass defined via fi = Xv 1 . The decay of the saxion into higgsinos 
is also efficient if it is allowed kinematically and its decay rate is given by 



1 rr? a ( ji 



8ir F% V m 



(15) 



In the DFSZ axion model, the decay mode into the MSSM particles may dominate over 
the two axion decay. The saxion decay temperature is then calculated as 

Note also that the mixing of the PQ scalar with Higgs fields induce the saxion decay into 
SM quarks and leptons. They are subdominant as long as the saxion is heavy (m CT > 
ITeV). 



3 Smooth hybrid inflation in a SUSY axion model 
3.1 The potential of the inflaton 

In this section, we review the smooth hybrid inflation model proposed originally in [12], 
including the supergravity correction. We start with the PQ sector superpotential (2) 
and the following Kahler potential, 

w = s (-^+j§S^)+ w °> ( i7 ) 

K= \S\ 2 + |^| 2 + |^| 2 + ..., (18) 

where ... in the Kahler potential denotes the higher order Planck suppressed terms. We 
have added a constant term Wo in the superpotential, which is required to cancel the 
positive vacuum energy from the SUSY breaking and make the cosmological constant 
zero in the present vacuum. The gravitino mass is related to it as Wq = m 3 / 2 Mp. We 
take Mp = 1 in this section. 

1 In our model, A must be an extremely small value such as A ~ 10~ 12 for v ~ 10 15 GeV in or- 
der to realize the appropriate /i-term. To avoid such a small A, we can adopt the superpotential 
W = \*i> m /M™~ 1 H d H d with some integer m > 2 instead of the second term in (2). Adopting such 
a superpotential, the decay rate of the saxion into the Higgs boson pair and that into the higgsino pair 
are m 2 times larger than (14) and (15). 



5 



The F-term scalar potential is calculated by using the formula 

V F = e K/M p [K ij *DiWD r W* - 3\W\ 2 /M 2 P ] , 



(19) 



where DjW = Wi + KiW/Mp and the subscript represents the derivative with respect to 
corresponding field and K % i* is the inverse matrix of Kij*. Then, the scalar potential is 
calculated as 



y =e is| 2 +l*l 2 +l*l 2 



;i + |5| 5 
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(20) 



-3 



Writing down only a few significant terms, the scalar potential is expressed as 



V =[ l+|^| 2 +| 1 |r| 2 + -| 1 9| 4 



A#2(n-1) 



+ 



n^Pl^p^-D „., ID|2 / (**) 



, r . , (|*| + 1*1) -2n/i 2 |S|" . 



+ c.c. 



(21) 



+ 2IU 



5 // 2 + (n-i; 



(^*) n 



+ c.c. 



+ . . . . 



If \S\ ^> (/iM n 1 ) 1 / n = / a is hold, the PQ fields are placed on the following temporal 
minimum 

f / f \ V(n-i) 

J a I J a 



1*1 = 1*1 



for 151 > f a . 



(22) 



[n(2n- l)]V2(n-i) 

Expressing the complex inflaton field by = \/2|5| and #s = arg(5) and substituting 
(22) into (21), we get the effective potential of the inflaton (p as 



Vifp) =H 4 



1 - 



4 (2n-l)/(n-l) / n _ 1 \/f a 



n — 1 



2n/(ra-l) 



+ 2y/2jj 2 m 3 / 2 (pcos6< : 



1 + 



fa 

[n{n - l/2)] n /( n - x ) 



+ > + ■ 



2n/(n-l)' 



(23) 



for |,S| 3> / . Under this potential, the inflation takes place along the temporal valley 
given by (22) for large field value of \S\. We numerically calculated the time evolutions of 
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Figure 1: The time evolution of Re(S) (red solid line), (green dashed line) and \^\ 
(blue dotted line) in the inflationary epoch are shown. The x-axis is the time multiplied 
by the Hubble parameter during inflation. We have taken [i = 0.002, M = 1.5 (f a = 0.05), 
n = 2 and W = and S = 0.3, * = 0.01 + 0.002z, ^ = 0.02 + O.OOlz and H = 2 x 10" 6 
as initial values in Planck units [Mp = 1). 



the fields Re(S'), and during and soon after the inflationary era, which are shown 
in Fig. 1. This supports the above analytical calculations. 

From the inflaton potential (23), the slow- roll parameters [20] are calculated as 



1 (V 
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(3n-l)/(n-l) 



(24) 



(25) 



where the prime denotes the derivative with respect to (p and (pa is defined as 

/ A 2n-l \ V(6n-4) 

<0 . = [ - I W(3n-2) 

The slow-roll condition is broken when \rj\ ~ 1 at which p reaches the value given by 



(26) 



n(3n 



(n - l)(2n - 1) 
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(27) 
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Figure 2: /x and function of f a are shown. The WMAP normalization of the 

density perturbation is imposed. We have taken n = 2 (Fig. 2(a) and 2(b)), n = 4 
(Fig. 2(c) and 2(d)), 9 S — 0, Wq — (solid-red line), m 3 / 2 = 10 TeV (dashed-green line), 
m 3/2 = 100 TeV (dotted-blue line) and m 3 / 2 = 1000 TeV (small-dotted-magenta line). 
The e-folding number is set to be 50 in all figures. 



Using these quantities, we can calculate the power spectrum of the curvature perturbation 
and scalar spectral index n s . Results are shown in Fig. 2. On each contour, the WMAP 
normalization of the density perturbation [5] is imposed. From these figures, the PQ 
breaking scale can be reduced to ~ 4 x 10 14 GeV and the spectral index well-agrees with 
the observational value. Fig. 3 shows the Hubble scale during inflation, Hj, as a function 
of f a by imposing the WMAP normalization. 
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(a) n=2 (b) n=4 

Figure 3: Hi as a function of f a after the WMAP normalization of the density pertur- 
bation is imposed. We have taken n = 2 in Fig 3(a) and n = 4 in Fig. 3(b) and Wq = 
(solid-red lines), m 3 / 2 = 10 TeV (dashed-green lines), m 3 / 2 = 100 TeV (dotted-blue lines) 
and 777,3/2 = 1000 TeV (small-dotted-magenta lines) and the e-folding number is set to be 
50 in both figures. 



3.2 The initial value of the inflaton 



In some parameter regions, the linear term in the potential of the inflaton (23) may create 
a local minimum if cos 9 s < 0, which may spoil the success of the inflation. 2 In such a 
case, the initial values of ip and 9s are constrained so that the inflaton is not trapped at 
the local minimum. Let us investigate these issues. 

First we examine whether a local minimum (and maximum) exists or not. Clearly there 
is no such minimum for cos 9 s > and hence we consider cos 9s < in the following. For 
notational simplicity, let us rewrite the inflaton potential (23) as 



V(<p) = V - A<p~ a -B<p + CV 
where Vo, A, B, C, a and (3 are defined as 

4 (2n-l)/(n-l) n _ 1 
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The derivative of the potential with respect to (p is given by 

-a-l / \ p-i 



V(<p) = B 



- 1 



(28) 



(29) 



(30) 



For the case of hybrid inflation, the effect of linear term was studied in Refs. [21, 22, 23] 
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(a) n = 2 (b) n = 4 

Figure 4: A parameters for m 3 / 2 = 10 TeV (solid red lines), 100 TeV (dashed green 
lines) and 1000 TeV (dotted blue lines) are shown. We have taken n = 2 (Fig. 4(a)) and 
n = 4 (Fig. 4(b)) and imposed the CMB normalization condition. 



where y9 max and (f m i n are defined as 

faA\ 1/{a+1) /fl^M / 

¥W = I "g- I , ¥?min = I ^ I , (31) 

corresponding to the local maximum for small C limit and the local minimum for small 
A limit, respectively. The condition for non-existing of the local maximum and the local 
minimum in the potential (28) is 

V'(<pv) > 0, (32) 
where ipy is defined through V"(ipy) = and given by 

„ = (<£±m v ™ ,33) 



m - i)c 



Thus the condition (32) is written as 



A /^^//^iy^w >L 

If this condition is satisfied, the local minimum of the potential does not appear and the 
initial value of 9s is not constrained. Moreover, there is no upper bound on the initial 
value of (p, and hence the initial value of the inflaton is only bound below so that inflation 
continues at least 50 e-foldings. We have calculated A under the various parameter sets 
which satisfy the CMB normalization and the results are shown in Fig. 4. We found that 
there is no local minimum for m 3 / 2 ^ 100 TeV (n = 2) or m 3 / 2 ^ 1000 TeV (n = 4) for 
/ a ~10 15 GeV. 
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From the above discussion, the local minimum of the inflaton potential may appear 
for large gravitino mass and relatively small n, which leads to A < 1. In such a case, the 
initial values of ip and 9s are constrained in order to avoid trapping at the local minimum. 
For successful inflation, one of the following conditions must be satisfied. First, the initial 
value of 9s must be small and the subsequent evolution must not change 9 s significantly 
in order to keep cos^s > where local minimum of ip does not appear. Second, the initial 
value of (p must be smaller than the local maximum given by 



2(5n-l)/(6n-2) j 
[2n{2n - l)]n/(3n-l) 



_(2n - l)m 3 / 2 /a(- cos 9 S ) 



(n-l)/(3n-l) 



(35) 



provided cos#g < 0. To investigate the allowed region of the initial value of the inflaton, 
we integrate the following equations of motion of the inflaton : 



3H(p ~ — /1 4 



4 (2n-l)/(n-l) 



n 



[2n{2n - 2n-l\(p J 



^yn/(n-l) 1 L 



<p 4 



3H9, 



2^2 ., 
A' 



fi 2 m 3/2 sin 6s 



m 3 / 2 cos9s, (36) 



(37) 



As a result, we show the allowed regions of the initial values of the inflaton field in 
Fig. 5. We have taken f a = 10 15 GeV, fi = 3 x 10 14 GeV and n = 2. From these results, 
there is no severe initial value problem in the smooth hybrid inflation model even if we 
adopt relatively large gravitino mass 777,3/2 ~ 100 TeV, which is contrasted to the case of 
SUSY hybrid inflation [23]. 

Before ending this section, we comment on the advantages of the smooth hybrid in- 
flation model compared with the hybrid inflation model [7, 8]. 

• Because the PQ symmetry is already broken during inflation, the initial misalign- 
ment angle at each spatial point takes an identical value at least in our observable 
Universe, and hence cosmic strings and domain walls are never formed after infla- 
tion. In the hybrid inflation model, on the other hand, the PQ symmetry is broken 
at the end of inflation and the initial misalignment angles take random values at 
different spacial points. In such a case, the cosmic strings are formed at the PQ sym- 
metry breaking and the domain walls are formed at the QCD phase transition [24], 
which induces the additional contribution to the axion cold dark matter [25, 26]. 

• In the smooth hybrid inflation model, an initial misalignment angle at each spatial 
point can be chosen to be arbitrarily small value from the same reason as above. 
This implies the PQ symmetry breaking scale can be larger than the usual upper 
limit 10 12 GeV without conflicting with the CDM density parameter even without 
late-time entropy production. 

• The scalar spectral index can be well consistent with the WMAP central value 
without invoking tuning in the non-minimal Kahler potential. In the hybrid inflation 
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(a) n = 2 and m 3/2 = 100 TeV (b) n = 2 and m 3/2 = 1000 TeV 

Figure 5: The allowed regions of the initial value of the inflaton are shown as the red 
regions. We have taken f a = 10 15 GeV, \i — 3 x 10 14 GeV and n = 2. The lower 
limit of if comes from the condition that the inflation continues at least 50 e-foldings. In 
Fig. 5(b), the dashed green line represennts Wax and the empty region corresponds to 
the breakdown of our assumption 9$ ~ constant. 



model, on the other hand, the scalar spectral index is relatively large (n s > 0.98), 
under the minimal Kahler potential. Although the WMAP central value can be 
reproduced under the non-minimal Kahler potential [27, 28], a severe initial value 
problem of the inflaton arises. These problems do not exist in the smooth-hybrid 
inflation model. 

• The presence of the constant term in the superpotential is not so effective com- 
pared with the SUSY hybrid inflation case. In the SUSY hybrid inflation model, 
the parameter space is strongly constrained for the large gravitino mass under the 
minimal Kahler potential [23]. In the smooth hybrid inflation model, we found 
that, from Fig. 2(a) and 2(b), the effects of the constant term is small up to m 3 / 2 ~ 
100 — 1000 TeV for most parameter regions. Note also that the angular motion of the 
inflaton field does not spoil the inflaton dynamics for 7713/2 < 100 — 1000 TeV even 
if 9s — vr/2 because the linear term in the inflaton potential can never dominate 
the dynamics, hence there is no severe initial value problem of the inflaton in the 
smooth hybrid inflation model for m 3 / 2 < 100 TeV. 

• A drawback of the smooth hybrid inflation model is that it predicts the large axion 
isocurvature perturbation, which poses severe constraint on the case of n = 2, as 
shown in Sec. 5. 
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Figure 6: k as a function of f a after the WMAP normalization of the density perturbation 
is imposed. We have taken n = 2 in Fig, 6(a) and n = 4 in Fig. 6(b) and Wo = (solid- 
red lines), m 3 / 2 = 10 TeV (dashed-green lines), m 3 / 2 = 100 TeV (dotted-blue lines) and 
m 3/2 = 1000 TeV (small-dotted-magenta line) and the e-folding number is set to be 50 in 
both figures. 

4 The scalar dynamics after inflation and cosmology 
4.1 Reheating after inflation 

In this section, we consider the PQ field dynamics after inflation. The post-inflationary 
dynamics of the PQ fields is determined by the following low-energy potential: 

2 
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+ 2m 



3/2 



jlf2(n-l) 

s(fi 2 + (n 
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n 2|£|2|^|2(n-l) 



M 4 (" _1 ) 



•(l*| 2 + 1*1 
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(38) 



+ ™3/ 2 (|S| 2 + l*| 2 + l*| 2 -3) 



Right after the inflation, the low-energy SUSY breaking terms written in the second row 



in (38) are negligible and the minima of the PQ fields are given by v 



f a - The flat 



direction which we identify as the saxion (a) and the direction perpendicular to the flat 
direction which is denoted by $ are related to the PQ scalar fields as 



1^1 



v 1 + 



a 



V2F a 



1^1 



v 1 + 



a 



V2F a 



(39) 



From the above relation, the first term in the right-hand-side of (38) gives the mass of $ 
and the second term gives the mass of S. $ and S have the same mass : m$ = mj = 
v2k/o for v = v = f a where k = n{f a /M) 2 ^ n ~ 1 ' which takes a value shown in Fig. 6. 
Furthermore they have the same initial amplitude of order f a , so they equally contribute 
to the total energy density of the universe. Because of the large mass and large coupling 
constant of $ and S, their lifetimes are quite short. 
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We must note the mixing of S and $. Redefining ip as Re(S)/\/2, the mixing term 
is written as 4n(n — l)m^/2fi 2 (f^/F a ~ m^^ms^ in the potential (38). Because of 
ms = m$, the mass eigenstates are (up ± $)/v / 2 whose squred masses are given by 
m 2 s ± 5m 2 respectively, where 5m 2 ~ m^ms- Hence the mixing time scale is estimated 
as 2ms /5m 2 ~ V m 3/2- This is much longer than the lifetime of both S and $ which is 
calculated below, so we can neglect the mixing effect here. 

Let us first consider the decay of S. From the superpotential (2), S has a Yukawa 
interaction with the axino. Representing the fermionic components of the PQ superfields 
\17 (^) as if) (if)) respectively, the axino field is defined as a = (if) — $)/v2, whose mass is 
given by /c|us|(= m 3 / 2 /n), where v$ = (S) = — 2m 3 / 2 /i 2 /m| = —m 3 / 2 /(nn) is the VEV of 
S. The other combinations of S (fermionic components of S) and (if> + if>)/ a/2 get masses 
of ~ y/2nf a . The interaction of the axino with S is derived from the superpotential (2) 
as 

— Cs-a = -uSaa + h.c. (40) 
From this term, the decay rate of S into axino pair is derived as 



V2 



Fs-^aa — QO — W^-^fa- (41) 
OZTT OZ7T 

In the present model, the coupling constant is relatively large such as K ~ 0.1 - 1, which 
implies the axino production occurs soon after inflation. The decay temperature of S 
defined through the axino energy density is calculated as 

/ 2 \ -!/ 4 / \ 3/2 / n \ 1/2 

r. (ts) =(^) ^S;.3xlO»Gev(Aj (42 ) 

where gs, = 2 x 7/8 = 1.75 is the relativistic degrees of freedom of axinos and ts represents 
the time of the S decay. Since the axino decoupling temperature is given by ~ 
10 17 GeV(/ a /10 15 GeV) 2 [29], the produced axinos cannot have any thermal contacts with 
the MSSM particles. For convenience, however, we call T s as the decay temperature of S. 

Now we consider the decay of $. Similar to the saxion, $ decays mainly into gluons 
and gluinos through X and X loops, so the Universe is reheated by the $ decay. The 
decay rate of $ into gluons is calculated as 

_^ 2 _m|_^ 2 _ 3 
^"" 327r3F a 2 "16^ /a ' [ 6) 

Thus the reheating temperature is given by 

2 \ - 1 / 4 / \ / \ 3/2 / n \ 1/2 



R 



hl g \ JT^ qq M P ~ 2 x 10 12 Gevf -^-) ( —) ( ^ 

V9o y V v 99 Vo.04 y Vo.i y Vio l5 Gevy 



, (44) 



where is the relativistic degree of freedom which we take g* ~ 100 in the last equality. 
After the decay of S into axinos until the decay of $, the universe is dominated by $ 
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because the produced axino from S is highly relativistic and its energy density is redshifted 
relative to that of $ as p„ oc a _1 p$. Thus we can calculate the axino abundance as 

= ~ 01 P° 15GeV l ( Tr ] (45) 

2m,- \k)\ f a Al0 12 GeV'' 



S 



which implies the overproduction of axinos, but the late-time entropy production due to 
the saxion decay dilutes it. Cosmological implications will be discussed below. 



4.2 Saxion dynamics and its implications 

After the decay of $, the Universe is radiation dominated. The existence of thermal 
background induces thermal effects on the scalar potential. In particular, in the KSVZ 
model, \I/ interacts with heavy quarks Q and Q and heavy quarks interact with with 
MSSM particles through the QCD couplings. Then the two-loop correction to the strong 
coupling constant induces the thermal effective potential [30] such as 

V th ^a 2 s THog^. (46) 

The saxion dynamics is similar to the previous studies [7, 8], and here we briefly repeat the 
discussion. This thermal- log potential lifts the flat direction, and (\^\) tends to roll 
down to the smaller (larger) value. Just after the reheating, both \1/ and ^ are placed at f a 
and the thermal mass is larger than the Hubble parameter. Thus rolls down towards 
a larger value without feeling the Hubble friction. The thermal mass decreases as \l/ rolls 
down and becomes comparable to the Hubble parameter and then \l/ gets frozen due to 
the Hubble overdamping at the field value of \^\ ~ a s Mp. As the Universe expands and 
temperature decreases, the thermal effect becomes insignificant and the zero-temperature 
minimum (7) appears. At this epoch, the PQ field is displaced from the zero-temperature 
minimum. As a result, when the Hubble parameter becomes comparable to the gravitino 
mass, the PQ field starts to oscillate around its true minimum along the flat direction, 
with an amplitude of a,i ~ a s Mp. 

In Fig. 7 we show the above mentioned dynamics by numerically solving the equation 
of motions of fields. In this figure, the time evolution of effective mass m t h{— a s T 2 /\^/\) 
(red solid line), Hubble parameter (green dashed line) and saxion field value (blue dotted 
line) are shown. The saxion energy density divided by the entropy density after the 
oscillation is thus given by 

Note that the entropy density in this expression includes only that produced by the 
inflaton decay. Since the saxion dominates the Universe before it decays, its decay dilutes 
the abundance of preexisting matter. The dilution factor is given by 

7 s ?Zk (*r' „ 2 x 10 -» (^-) ( HlfXy 72 (48) 

1 4 \sJ VlOOMeVy V rn a J \ ai J y ' 
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Figure 7: The dynamics of the saxion after reheating are shown. We show the time 
evolution of m t h (red solid line), Hubble parameter (green dashed line) and saxion (blue 
dotted line). We have taken /i = 0.05, M = 0.1 (f a = 0.07)_, a s = 1, m 3/2 = 10~ 5 , c x = 1 
and C2 = 2 in units of Mp = 1 and S = 0.1f a and \& = \I/ = / a as initial values. The 
small-dotted magenta line represents the gravitino mass. 



Number densities per entropy density of all matters, which are produced at the reheating 
and conserved thereafter such as the gravitino and axino, are diluted by this factor. As a 
result, thermally produced gravitinos [31] as well as the non-thermally produced ones by 
the direct inflaton decay [32, 33] are also diluted away to a cosmologically negligible level. 
Thermally produced axino abundance during the reheating [34, 35] is also negligibly small 
due to the large PQ scale and huge dilution. 

On the other hand, the axino abundance produced by the inflaton decay, given by 
Eq. (45), may not be neglected even after the dilution by the saxion decay. The axino 
obtains a mass of ~ iTL 3 / 2 / n and its decay width into the gluino and gluon is given by 

^Lul (A9 \ 

^gg- 167r 3 F 2- V*) 

Since the decay width is comparable to the saxion in the KSVZ model, the LSPs produced 
by the axino decay is always smaller than or comparable to those produced from the saxion 
decay. 3 Below we examine the LSP abundance from the saxion decay, taking account of 
self-annihilation of the LSP. 



3 In the DFSZ model, the axino decay rate into Higgs boson and higgsino is comparable to the saxion 
decay rate into the higgsino pair (15). 
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4.3 Wino dark matter from saxion decay 

In this section we see that Winos produced by the saxion decay can account for the present 
dark matter abundance. 4 Assuming the anomaly mediated SUSY breaking (AMSB) 
model [36, 37], or the pure-gravity mediation model [38] in which the 125 GeV Higgs 
boson can easily be explained, the gravitino mass is of the order of ~ 10 2 - 10 3 TeV. The 
AMSB contributions to the the MSSM gaugino masses, M a = (Mi, M 2 , M 3 ), are given by 

Ma = lfe^ m3/2 ' (50) 

where g a are the SM gauge coupling constants and b a = (11, 1, —3). Then, the ratio of 
the MSSM gaugino masses are given by : : m- g ~ 3 : 1 : 8 implying the Wino- 
LSP with mass of 0(100) GeV - 0(1) TeV, although the relation may be modified by the 
Higgs-higgsino loop contribution. Because of the large gravitino mass, the saxion decay 
temperature given by (13) or (16) becomes 0(1) GeV. In such a case, the present dark 
matter abundance can be explained by the Wino-LSP from the saxion decay as shown 
below. 

Because we consider the large gravitino mass of order 0(100) TeV and the LSP mass 
is given by m x ~ m 3/ / 2 /400, the freeze-out temperature of the LSP, given by T fr « m x /25, 
may be higher than the saxion decay temperature, which implies the produced LSPs 
are never in thermal equilibrium. Hence, in the following arguments, we focus on the 
non-thermally produced dark matter [39, 40, 19]. 

Since the decay rate of the saxion into SUSY particles is comparable to that into 
ordinary particles, (see (11) and (12) in the KSVZ model or (14) and (15) in the DFSZ 
model) a large number of LSPs are produced. The annihilation cross section of Wino-LSP 
is given by [39] 

, gi 1 (l-x K )W 

M = 2^f (2 -*»,)» ■ (51) 

where Xw = mw/ m x with mw being the W boson mass. This is roughly estimated as 
(vannv) ~ 10~ 7 GeV~ 2 for m x ~ 100 GeV. Note that the Wino cross section is constrained 
from WMAP observations [41, 42] and recent gamma-ray measurements by the Fermi 
satellite [43], which excludes the Wino mass below ~ 300 GeV as a dominant component 
dark matter. Note that we have ignored the Sommerfeld enhancement effect [44] . This is 
valid as long as we consider the Wino lighter than ~ 1 TeV [45] . 

The abundance of the LSP is determined by solving the following Boltzmann equation: 

diTh 

+ 3Hn x = -(a ann v)n 2 x , (52) 
where n x is the LSP number density. Assuming (<J ann v) is temperature-independent, this 



4 In the next section, we will see that the axion cannot be the dominant component of DM once the 
isocurvature constraint is taken into account. 
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Figure 8: The relations between the wino abundance and the saxion mass are shown. The 
solid red lines, the dashed green lines and dotted blue line correspond to the KSVZ model, 
the DFSZ model with /i = 50m x (higgsino mass) and the DFSZ model with fi = 100m x 
respectively. We have taken m a = ^3/2 (thick lines) and m a = 2m^/2 (thin lines). The 
breakiing points reflect the constraint on the Wino mass, m x > 300 GeV. 



is analytically solved and we get the LSP abundance in terms of Y X (T) = n x (T)/s as 

/o 2 \ 1/2 

Y-\T) = Y~\T a ) + [j^9*{T„)) (a ann v)M P (T a - T), (53) 

where Y x (T a ) contains the contributions from both non-thermally produced LSPs by the 
saxion and axino decay and also thermal freeze-out LSPs. From this expression, we 
soon realize that the result is independent of the initial abundance Y x (T a ) if the initial 
abundance or the annihilation cross section is large enough, as is expected. Actually, 
in the present setup, the abundance of LSPs from the saxion decay is large enough to 
annihilate efficiently. Thus the resultant density parameter is estimated as 

ni2 / 60 \ 1/2 / m x \ /KT 8 GeV- 2 \ ZlGeV\ , . 

This implies that the non-thermally produced Winos can eventually becomes the dominant 
component of the dark matter. Taking into account the cosmological and astrophysical 
constraint on the Wino mass, m x > 300 GeV [38], the resultant abundance of the Wino 
dark matter is shown in Fig. 8. From this, the mass of the saxion is predicted as m a ~ 
100 TeV in both KSVZ and DFSZ axion models with fi > 50m x . Such a Wino LSP may 
be found by direct or indirect detection experiments or LHC [46, 47]. 
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5 Isocurvature perturbation and non-Gaussianity 



In this section, we calculate the axion abundance and the constraints from the isocurvature 
perturbation. We also calculate the non-Gaussianity from the isocurvature perturbation. 

We first note that the axion cannot be diluted by the saxion decay once we impose 
T CT > 1 GeV in order to realize the Wino dark matter as shown above. Hence we have to 
tune the initial misalignment angle to suppress the axion abundance. The present density 
parameter of the axion is given by [48, 49, 50] 



Q n h z 



0.2 
0.2 



fa 



10 12 GeV 

fa 

10 12 GeV 



1.18 



o 1 



-0.82 



Hjf a /27r\^(k )\ 
10 12 GeV 



for 9i > 
for Qi < 



Hi 



2tt|^(A;o)| 
2tt|^(A; )| ! 



(55) 



which means that the axion abundance can no longer be reduced by tuning 9i if 0, is 
smaller than the critical value given by 



9 r 



# 7 /27T 



2 x 10" 4 



#//|*(*o)| 



10 



-3 



(56) 



where ^f(ko) represents the PQ field value when the pivot scale leaves the horizon during 
inflation. 5 



5.1 Constraint from isocurvature perturbation 

Since the PQ symmetry is broken during inflation, there exists the quantum fluctuation 
of the axion, so the large CDM isocurvature perturbation can be induced from the axion 
and the model parameters are strongly constrained [51, 52]. The ratio of the adiabatic 
curvature perturbation to the isocurvature one is parameterized as 

where and Vs are dimensionless power spectrum of the curvature and CDM isocur- 
vature perturbations, and a is constrained from the observation as a < 0.077 [5]. The 
power spectrum of the CDM isocurvature perturbation is related to that of the axion 
isocurvature perturbation in terms of the density parameters as 

a x 2 



'CDM 



V s = V S ,a (58) 

J if 



5 Since H/^ decreases during inflation, the critical value 8 CI should be evaluated when the largest 
observable scale (~ kg) leaves the horizon. 
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Figure 9: The power spectrum of the isocurvature perturbation (Fig. 9(a)) and the non- 
linearity parameter (Fig. 9(b)) are shown. In both figures, we have taken m 3 / 2 = 100 TeV 
and n = 2 (solid red line), n = 3 (dashed green line) and n = 4 (dotted blue line). 
In the left figure, the dashed-and-dotted cyan line represents the upper limit from the 
observation. 



and the power spectrum of the axion isocurvature perturbation is calculated as 



S.a 



Hr 



ir\*(ko)\0i 

1 

4 



for 9i > 
for 9i < 



27r|*(A;o)r 

gz 
2n\^(k )\ 



(59) 



and 

Substituting (55) and (59) into (58) and using the WMAP best fit values V ( = 2.43 x 
10~ 9 and ficDM^ 2 — 0.112 [5], the power spectrum the CDM isocurvature perturbation is 
rewritten as 



V 



1/2 



20,; 



fa 

10 15 GeV 



1.18 



Hi/\*(ko)\ 



4 x 10- 



fa 



10 15 GeV 



io- 3 

1.18 / u l\i\<n,. M\ 2 



10 



-3 



for 6> > 



for 6>, < 



2tt|^(A;o)| 

Hr 
27r|^(fco)| 



(60) 



This implies obviously that Qi > if//27r|\l/(fco)| is not allowed by the constraint from 
the observation and the small initial misalignment angle satisfying \^(ko)\6i < Hj/2n is 
necessary for the case of f a ~ 10 15 GeV. As seen from Eq. (55), the axion cannot be 
the dominant component of the CDM. The constraint on the parameters are shown in 
Fig. 9(a), in which 9^ < if//27r|\I/(/co)| is assumed. It is found that n < 3 is forbidden by 
the observation but n = 4 is allowed for some range of f a of our interest. 
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5.2 Non-Gaussianity from isocurvature perturbation 

In the inflationary paradigm, the primordial density perturbation, which is originated 
from the quantum fluctuation of the scalar field, obeys the almost Gaussian statistics. 
Thus if the deviation from the Gaussianity is detected, it has rich information on the 
origin of primordial density perturbation. In particular, the axion model can generate 
sizable non- Gaussianity in the isocurvature perturbation [54, 53]. The non-Gaussianity 
is characterized by the non-linearity parameter /nl defined by 

(C(kMh)C(h)) = pM^f^ih + k 2 + k 3 )[Pd k i)Pdk2) + (2 perms)], (61) 

where '2 perms' means 2 permutations. ( denote the curvature perturbation evaluated 
on the uniform density slicing and (k) is the power spectrum of ( defined as 

(((k)((k')) = (2nfP ( (k)6^(k + k'). (62) 

In the present model, the isocurvature perturbation from the axion becomes large and 
the axion density parameter is forced to be extremely small due to the constraint from 
the isocurvature perturbation, so a relatively large non-Gaussianity is expected. In the 
case of \^(k Q )\9i < Hj/271, the non-linearity parameter in our model is calculated as [53] 

/r = I^p(^) 3/2[ln(fcL)rl/2 ' (63) 

where L is an infrared cutoff which is taken to be the present Hubble horizon scale and 
we set [h^/cL)] 1 / 2 = 5 [55]. The /nl - fa plane is shown in Fig. 9(b), from which a 
relatively large non-Gaussianity with /nl ~ 0.1 - 1 is realized. It should be noticed that 
the we presented here is the non-Gaussianity in the (uncorrelated) CDM isocurvature 
perturbation [53, 56, 57, 58, 59], whose properties are different from non-Gaussianity in 
the adiabatic perturbation. 

6 Conclusions 

We have shown that the smooth hybrid inflation naturally takes place in a SUSY axion 
model in which the PQ fields are identified with a part of the inflaton sector. In order to 
reproduce the WMAP observation of the density perturbation, the PQ symmetry breaking 
scale must be of order of 10 15 GeV. The spectral index can naturally be the WMAP best 
fit value within the minimal Kahler potential. Because the PQ symmetry is already 
broken during inflation, topological defects such as the cosmic strings and domain walls 
are never formed in this model. After the inflation, the Universe is reheated by the decay 
of the heavy fields into the ordinary particles. We have followed the dynamics of light 
scalar fields, saxion, after inflation, and found that the saxion starts to oscillate with large 
initial amplitude of order a s Mp. Thus the saxion eventually dominates the Universe and 
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the decay of the saxion produces huge amount of entropy, which dilutes the gravitinos 
and axinos produced during reheating. It can also produce Wino DM nonthermally with 
correct dark matter abundance if the Wino is much lighter than the gravitino, as is 
expected from the AMSB or pure-gravity mediation model. A severe constraint is imposed 
in this model due to the axion isocurvature perturbation, since the PQ symmetry is 
broken during inflation. It excludes the possibility of axion coherent oscillation as a 
dominant component of the current CDM. A non-Gaussianity of the isocurvature type, 
/nl°' ) ~ 0.1 — 1, is predicted in our model, which may be detected by the future observation. 

Acknowledgment 

This work is supported by Grant-in- Aid for Scientific research from the Ministry of Edu- 
cation, Science, Sports, and Culture (MEXT), Japan, No. 14102004 (M.K.), No. 21111006 
(M.K. and K.N.), No. 22244030 (K.N.) and also by World Premier International Research 
Center Initiative (WPI Initiative), MEXT, Japan. N.K. is supported by the Japan Society 
for the Promotion of Science (JSPS). 

References 

[1] G. Aad et al. [ATLAS Collaboration], Phys. Lett. B 716, 1 (2012) [arXiv: 1207.7214 
[hep-ex]]; S. Chatrchyan et al. [CMS Collaboration], Phys. Lett. B 716, 30 (2012) 
[arXiv:1207.7235 [hep-ex]]. 

[2] C. A. Baker, D. D. Doyle, P. Geltenbort, K. Green, M. G. D. van der Grinten, 
P. G. Harris, P. Iaydjiev and S. N. Ivanov et al, Phys. Rev. Lett. 97, 131801 (2006) 
[hep-ex/0602020]. 

[3] R. D. Peccei and H. R. Quinn, Phys. Rev. Lett. 38, 1440 (1977). 

[4] For reviews, see J. E. Kim, Phys. Rept. 150, 1 (1987); J. E. Kim and G. Carosi, Rev. 
Mod. Phys. 82, 557 (2010). 

[5] E. Komatsu et al. [WMAP Collaboration], Astrophys. J. Suppl. 192, 18 (2011) 
[arXiv:1001.4538 [astro-ph.CO]]. 

[6] For a review, see S. P. Martin, arXiv:hep-ph/9709356. 

[7] M. Kawasaki, N. Kitajima and K. Nakayama, Phys. Rev. D 82, 123531 (2010) 
[arXiv: 1008.5013 [hep-ph]]. 

[8] M. Kawasaki, N. Kitajima and K. Nakayama, Phys. Rev. D 83, 123521 (2011) 
[arXiv: 1104. 1262 [hep-ph]]. 



22 



[9] E. J. Copeland, A. R. Liddle, D. H. Lyth, E. D. Stewart and D. Wands, Phys. Rev. D 
49, 6410 (1994) [astro-ph/9401011]; G. R. Dvali, Q. Shall and R. K. Schaefer, Phys. 
Rev. Lett. 73, 1886 (1994) [hep-ph/9406319]; A. D. Linde and A. Riotto, Phys. Rev. 
D 56, 1841 (1997) [hep-ph/9703209]. 

[10] P. J. Steinhardt and M. S. Turner, Phys. Lett. B 129, 51 (1983); G. Lazarides, 
C. Panagiotakopoulos and Q. Shan, Phys. Lett. B 192, 323 (1987); G. Lazarides, 
R. K. Schaefer, D. Seckel and Q. Shah, Nucl. Phys. B 346, 193 (1990). 

[11] M. Kawasaki, T. Moroi and T. Yanagida, Phys. Lett. B 383, 313 (1996) [hep- 
ph/9510461]. 

[12] G. Lazarides and C. Panagiotakopoulos, Phys. Rev. D 52, 559 (1995) [hep- 
ph/9506325]. 

[13] M. Yamaguchi and J. 'i. Yokoyama, Phys. Rev. D 70, 023513 (2004) [hep- 
ph/0402282]. 

[14] J. E. Kim, Phys. Rev. Lett. 43, 103 (1979); M. A. Shifman, A. I. Vainshtein and 
V. I. Zakharov, Nucl. Phys. B 166, 493 (1980). 

[15] M. Dine, W. Fischler and M. Srednicki, Phys. Lett. B 104, 199 (1981); A. R. Zhit- 
nitsky, Sov. J. Nucl. Phys. 31, 260 (1980) [Yad. Fiz. 31, 497 (1980)]. 

[16] T. Kugo, I. Ojima, T. Yanagida, Phys. Lett. B135, 402 (1984). 

[17] E. J. Chun and A. Lukas, Phys. Lett. B 357, 43 (1995) [hep-ph/9503233]. 

[18] M. Kawasaki, K. Nakayama and M. Senami, JCAP 0803, 009 (2008) 
[arXiv:0711.3083 [hep-ph]]. 

[19] M. Endo and F. Takahashi, Phys. Rev. D 74, 063502 (2006) [hep-ph/0606075]. 

[20] A. R. Liddle and D. H. Lyth, "Cosmological Inflation and Large-scale Structure," 
(Cambridge University Press, Cambridge, England, 2000). 

[21] W. Buchmuller, L. Covi and D. Delepine, Phys. Lett. B 491, 183 (2000) [hep- 
ph/0006168]. 

[22] V. N. Senoguz and Q. Shan, Phys. Rev. D 71, 043514 (2005) [hep-ph/0412102]. 

[23] K. Nakayama, F. Takahashi and T. T. Yanagida, JCAP 1012, 010 (2010) 
[arXiv: 1007.5152 [hep-ph]]. 

[24] P. Sikivie, Phys. Rev. Lett. 48, 1156 (1982). 

[25] T. Hiramatsu, M. Kawasaki, T. Sekiguchi, M. Yamaguchi and J. 'i. Yokoyama, Phys. 
Rev. D 83, 123531 (2011) [arXiv:1012.5502 [hep-ph]]. 



23 



[26] T. Hiramatsu, M. Kawasaki, K. 'i. Saikawa and T. Sekiguchi, Phys. Rev. D 85, 105020 
(2012) [Erratum-ibid. D 86, 089902 (2012)] [arXiv:1202.5851 [hep-ph]]; T. Hiramatsu, 
M. Kawasaki, K. 'i. Saikawa and T. Sekiguchi, arXiv: 1207.3166 [hep-ph]. 

[27] M. Bastero-Gil, S. F. King and Q. Shafi, Phys. Lett. B 651, 345 (2007) [hep- 
ph/0604198]. 

[28] M. ur Rehman, V. N. Senoguz and Q. Shafi, Phys. Rev. D 75, 043522 (2007) [hep- 
ph/0612023]. 

[29] K. Rajagopal, M. S. Turner and F. Wilczek, Nucl. Phys. B 358, 447 (1991). 

[30] A. Anisimov and M. Dine, Nucl. Phys. B 619, 729 (2001) [hep-ph/0008058]. 

[31] M. Bolz, A. Brandenburg and W. Buchmuller, Nucl. Phys. B 606, 518 (2001) 
[Erratum-ibid. B 790, 336 (2008)] [hep-ph/0012052]; J. Pradler and F. D. Steffen, 
Phys. Rev. D 75, 023509 (2007) [hep-ph/0608344]; Phys. Lett. B 648, 224 (2007) 
[hep-ph/0612291]; V. S. Rychkov and A. Strumia, Phys. Rev. D 75, 075011 (2007) 
[hep-ph/0701104]. 

[32] M. Endo, K. Hamaguchi and F. Takahashi, Phys. Rev. Lett. 96, 211301 (2006) 
[hep-ph/0602061]; S. Nakamura and M. Yamaguchi, Phys. Lett. B 638, 389 (2006) 
[hep-ph/0602081]; M. Dine, R. Kitano, A. Morisse and Y. Shirman, Phys. Rev. D 73, 
123518 (2006) [hep-ph/0604140]; M. Endo, K. Hamaguchi and F. Takahashi, Phys. 
Rev. D 74, 023531 (2006) [hep-ph/0605091]. 

[33] M. Kawasaki, F. Takahashi and T. T. Yanagida, Phys. Lett. B 638, 8 (2006) 
[hep-ph/0603265]; Phys. Rev. D 74, 043519 (2006) [hep-ph/0605297]; T. Asaka, 
S. Nakamura and M. Yamaguchi, Phys. Rev. D 74, 023520 (2006) [hep-ph/0604132]; 
M. Endo, F. Takahashi and T. T. Yanagida, Phys. Lett. B 658, 236 (2008) 
[hep-ph/0701042]; Phys. Rev. D 76, 083509 (2007) [arXiv:0706.0986 [hep-ph]]; 
K. Nakayama, F. Takahashi and T. T. Yanagida, arXiv: 1209.2583 [hep-ph]. 

[34] L. Covi, H. -B. Kim, J. E. Kim and L. Roszkowski, JHEP 0105, 033 (2001) [hep- 
ph/0101009]; A. Brandenburg and F. D. Steffen, JCAP 0408, 008 (2004) [hep- 
ph/0405158]; A. Strumia, JHEP 1006, 036 (2010) [arXiv:1003.5847 [hep-ph]]. 

[35] E. J. Chun, Phys. Rev. D 84 (2011) 043509 [arXiv:1104.2219 [hep-ph]]; K. J. Bae, 
K. Choi and S. H. Im, JHEP 1108, 065 (2011) [arXiv: 1106.2452 [hep-ph]]; K. J. Bae, 
E. J. Chun and S. H. Im, JCAP 1203, 013 (2012) [arXivrl 111. 5962 [hep-ph]]. 

[36] L. Randall and R. Sundrum, Nucl. Phys. B 557, 79 (1999) [hep-th/9810155]. 

[37] G. F. Giudice, M. A. Luty, H. Murayama and R. Rattazzi, JHEP 9812, 027 (1998) 
[hep-ph/9810442]. 



24 



[38] M. Ibe and T. T. Yanagida, Phys. Lett. B 709, 374 (2012) [arXiv: 11 12.2462 [hep- 
ph]]; M. Ibe, S. Matsumoto and T. T. Yanagida, Phys. Rev. D 85, 095011 (2012) 
[arXiv: 1202.2253 [hep-ph]]. 

[39] T. Moroi and L. Randall, Nucl. Phys. B 570, 455 (2000) [hep-ph/9906527]. 

[40] M. Fujii and K. Hamaguchi, Phys. Lett. B 525, 143 (2002) [hep-ph/0 110072]; Phys. 
Rev. D 66, 083501 (2002) [hep-ph/0205044]. 

[41] T. R. Slatyer, N. Padmanabhan and D. P. Finkbeiner, Phys. Rev. D 80, 043526 (2009) 
[arXiv:0906.1197 [astro-ph.CO]]; D. P. Finkbeiner, S. Galli, T. Lin and T. R. Slatyer, 
Phys. Rev. D 85, 043522 (2012) [arXiv: 1109.6322 [astro-ph.CO]]. 

[42] T. Kanzaki, M. Kawasaki and K. Nakayama, Prog. Theor. Phys. 123, 853 (2010) 
[arXiv:0907.3985 [astro-ph.CO]]; J. Hisano, M. Kawasaki, K. Kohri, T. Moroi, 
K. Nakayama and T. Sekiguchi, Phys. Rev. D 83, 123511 (2011) [arXiv:1102.4658 
[hep-ph]]. 

[43] For recent analysis, see e.g., D. Hooper, C. Kelso and F. S. Queiroz, arXiv: 1209.3015 
[astro-ph.HE]. 

[44] J. Hisano, S. Matsumoto and M. M. Nojiri, Phys. Rev. Lett. 92, 031303 (2004) [hep- 
ph/0307216]; J. Hisano, S. .Matsumoto, M. M. Nojiri and O. Saito, Phys. Rev. D 
71, 063528 (2005) [hep-ph/0412403]. 

[45] J. Hisano, S. Matsumoto, M. Nagai, O. Saito and M. Senami, Phys. Lett. B 646, 34 
(2007) [hep-ph/0610249]. 

[46] T. Moroi and K. Nakayama, Phys. Lett. B 710, 159 (2012) [arXiv:1112.3123 [hep- 
ph]]. 

[47] B. Bhattacherjee, B. Feldstein, M. Ibe, S. Matsumoto and T. T. Yanagida, 
arXiv: 1207.5453 [hep-ph]. 

[48] D. J. Gross, R. D. Pisarski and L. G. Yaffe, Rev. Mod. Phys. 53, 43 (1981). 

[49] M. S. Turner, Phys. Rev. D 33, 889 (1986). 

[50] D. H. Lyth, Phys. Rev. D 45, 3394 (1992). 

[51] R. Bean, J. Dunkley and E. Pierpaoli, Phys. Rev. D 74, 063503 (2006) [astro- 
ph/0606685]. 

[52] M. Kawasaki and T. Sekiguchi, Prog. Theor. Phys. 120, 995 (2008) [arXiv:0705.2853 
[astro-ph]]. 

[53] M. Kawasaki, K. Nakayama, T. Sekiguchi, T. Suyama and F. Takahashi, JCAP 0811, 
019 (2008) [arXiv:0808.0009 [astro-ph]]. 



25 



[54] L. Boubekeur and D. .H. Lyth, Phys. Rev. D 73, 021301 (2006) [astro-ph/ 05 04046]. 

[55] D. H. Lyth, JCAP 0712, 016 (2007) [arXiv:0707.0361 [astro-ph]]. 

[56] M. Kawasaki, K. Nakayama, T. Sekiguchi, T. Suyama and F. Takahashi, JCAP 0901, 
042 (2009) [arXiv:0810.0208 [astro-ph]]; E. Kawakami, M. Kawasaki, K. Nakayama 
and F. Takahashi, JCAP 0909, 002 (2009) [arXiv:0905.1552 [astro-ph.CO]]; 
E. Kawakami, M. Kawasaki, K. Miyamoto, K. Nakayama and T. Sekiguchi, JCAP 
1207, 037 (2012) [arXiv: 1202.4890 [astro-ph.CO]]. 

[57] D. Langlois, F. Vernizzi and D. Wands, JCAP 0812, 004 (2008) [arXiv: 0809. 4646 
[astro-ph]]; D. Langlois and A. Lepidi, JCAP 1101, 008 (2011) [arXiv: 1007.5498 
[astro-ph.CO]]; D. Langlois and B. van Tent, Class. Quant. Grav. 28, 222001 (2011) 
[arXiv: 1104.2567 [astro-ph.CO]]; D. Langlois and B. van Tent, JCAP 1207, 040 
(2012) [arXiv: 1204.5042 [astro-ph.CO]]. 

[58] C. Hikage, K. Koyama, T. Matsubara, T. Takahashi and M. Yamaguchi, Mon. 
Not. Roy. Astron. Soc. 398, 2188 (2009) [arXiv:0812.3500 [astro-ph]]; C. Hikage, 
D. Munshi, A. Heavens and P. Coles, Mon. Not. Roy. Astron. Soc. 404, 1505 (2010) 
[arXiv:0907.0261 [astro-ph.CO]]. 

[59] C. Hikage, M. Kawasaki, T. Sekiguchi and T. Takahashi, arXiv:1211.1095 [astro- 
ph.CO]. 



26 



